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Abstract 

We analyze the atmospheric neutrino experiments of Super-Kamiokande (830- 
920 live days) using the three-flavor neutrino framework with the mass hierarchy 
mi ~ ni2 <C m^. We study the sub-GeV, multi-GeV neutrinos and upward through- 
going and stopping muons zenith angle distributions taking account of the Earth 
matter effects thoroughly. We obtain the allowed regions of mass and mixing pa- 
rameters Am^, #i3 and 623- In our present analysis, we used the solar neutrino 
small angle solution and large angle solution for Am 2 2 and 9\2- Am^ is restricted 
to 0.002-0.01 eV 2 and 9 13 < 13°, 35° < 6» 23 < 55° in 90% C.L. For 9 12 , there is 
no difference between the large angle solution and the small one. From x 2 fit; the 
minimum x 2 = 55 (54 DOF) is obtained at Ara^ = 4 x 10~ 3 eV 2 , 613 = 10° and 
#23 = 45°. In a two flavor mixing approximation (#13 = 0), the minimum x 2 = 61 
(54 DOF) is obtained at Am| 3 = 3 x 10~ 3 eV 2 and 6» 23 = 45°. If 13 = 10° is real, 
the detected v e events in K2K experiment will be about 10 times as large as events 
expected in #13 = case . 
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1 INTRODUCTION 



Super-Kamiokande collaboration has confirmed a neutrino oscillation by their atmospheric 
neutrino experiments In their two-flavor mixing analyses of the sub-GeV and the 
multi-GeV zenith angle distribution, it has been obtained that the <-» v T oscillation is 
preferred to the <-> v e oscillation and the range of mass parameter Am 2 is from 10 _3 eV 2 
to 10~ 2 eV 2 . Recently, Super-Kamiokande collaboration has reported that the <-> v T 
oscillation is maximally sin 2 26 = 0.9-1 and mass parameter Am 2 is from 2.5 x 10 _3 eV 2 
to 5 x 10 _3 eV 2 , using the sub-GeV, multi-GeV neutrino and upward muons zenith angle 
distribution experiments (830-920 live days) ||. 

However, these results are obtained from the two-flavor neutrino analyses. In order to 
account for the solar neutrino anomaly data together with atmospheric neutrino experi- 
ments, three flavor neutrinos are necessary at least. In three-flavor neutrinos scenario 0, 
with a mass hierarchy mi ~ m 2 <C m 3 , there are necessary two mass parameters Am 2 2 and 
Am 23 , and three mixing angles 6\2, 613 and 823- Solar neutrino anomaly gives constraint on 
only three parameters Am 2 2 , # 12 and 613. The MSW solution for solar neutrinos predicts 
the large mixing angle solution (Am 2 2 = 4 x 10~ 6 - 7 x 10~ 5 eV 2 , sin 2 29i 2 = 0.6-0.9) and 
the small mixing angle solution (Am 2 2 = 3 x 10~ 6 -1.2 x 10~ 5 eV 2 , sin 2 26 l2 = 0.003-0.01) 
for 9\3 = 0°-20°, and these large and small mixing angle solutions are merged for 25°-50° 
||. The vacuum solution is also obtained as Am 2 2 ~ 10 _10 eV 2 . However, CHOOZ experi- 
ment which is a terrestrial experiment using reactor neutrinos gives a strong constraint 
sin 2 2^13 < 0.18 for large Am 23 and Am 23 < 0.9 x 10~ 3 eV 2 for sin 2 26>i 3 ~ 1. From the 
recent many atmospheric neutrino analyses [[|, 0, [7], Q using either two- or three-flavor 
neutrinos framework, Am 23 = 10 _3 -10~ 2 eV 2 is obtained. Therefore, it can be said that 
the mixing angle 9 13 is small but not necessarily 0. 

In the three-flavor neutrino framework, if the mixing parameter 613 is not zero, it is 
necessary to consider the interplay between two mass parameters Am 2 2 , Am 23 and three 
mixing angles 612, #i 3 , 623 for atmospheric neutrino experimental analyses. The oscillation 
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term sin 2 1.27 — ^r 2 "-^ cannot be neglected though Am 2 2 <C Am| 3 , because the term is not 
so small in the sub-GeV experiment (E = 0.2 ~ 1.3GeV) of atmospheric neutrino with 
zenith angle 9 ~ 180° at which L ~ 10000km. In the multi-GeV experiment (E > 1.3GeV) 
of atmospheric neutrino and terrestrial short- and long-baseline experiments, the neutrino 
oscillation term sin 1. can be neglected. The angle #13 has been seen to be small, 

then the #13 is approximated to be in usual analyses. However, in the case considering 
the matter effects of the Earth, above approximation is not appropriate. The matter effect 
is expressed by the induced mass squared A = 2^/2EG F N e = 7.59pE x 10- 5 eV 2 , where 
N e is the number density of electrons, E is the energy of neutrinos (measured in GeV) 
and p is the matter density (measured in gm/cm 3 ). The Earth density p is 3.5 ~ 13g/cm 3 
and the multi-GeV experiment energy of neutrinos ranges from 1.3GeV to lOOGeV. Then 
the value of A go through the value of mass parameter Am^ and the MSW effect occurs 
in the mixing angle #13 (see section 3). 

In this paper, we analyze the Super-Kamiokande atmospheric neutrino experimental 
data(850-920 live days) of sub-GeV and multi-GeV neutrino, upward through-going 
and stopping muons, using the three-flavor neutrino framework with the hierarchy rrii m 
m 2 <C m 3 . In previous paper ||, we analyzed the 535 live days data of Super-Kamiokande 
atmospheric neutrino experiment |l[] using similar framework to the present analysis, 
approximating the Earth density to be constant. In this present analysis, we will pursue 
a full calculation adopting the varying density of the Earth. We will also discuss the 
long-baseline K2K experiment |§ in the three-flavor neutrino framework. 

2 NEUTRINO OSCILLATION IN THREE-FLAVOR 
NEUTRINOS 

The unitary matrix U which transforms the mass eigenstate neutrinos v a to the flavor 

eigenstate neutrinos z/; as the formula 

3 

Vi = U laVa, l = e, /I, r, (1) 
0=1 
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is parametrized as follows: 

U = exp (1623X7) exp (i9 13 X 5 ) exp (id 12 X 2 ) 

^ C12C13 S12C13 S13 

= — S12C23 — C 12 S 2 3Si3 C12C23 — S 12 S23Si3 S23 c 13 | , (2) 
\ S12S23 — C 12 C23Si3 — C12S23 — S 12C23 S 13 C23 c 13 , 

— cos $ ij 5 — sin $ %j • 

in a case disregarding the CP violation. We assume the mass hierarchy 

m\ ~ m 2 <C m 3 , (3) 

then Am| 2 <C Am 2 3 ~ Am| 3 . In this mass hierarchy, the transition probabilities P(v\ — > 
z///) can be written as 

P(n-n) = 1 - 2(1 - 2Ul - Ufi - U* + U? 3 )S 12 - 4Uf 3 (l - Uf 3 )S 23 , (4) 
P(u^u v ) = P(i^/ - i^)= 2(1^^ + ^^-^3)5^+ 41^^38, (5) 

where is a term representing the neutrino oscillation defined as; 

S a , = sin 2 1.27^L. (6) 

Here AmL = |m 2 — P and L are measured in units eV 2 , GeV and km, respectively. 
For the mass parameter Amf 2 and mixing angle #12, we use values obtained in the solar 
neutrino experimental analyses ||: the large mixing angle solution Am 2 2 = 3 x 10 _5 eV 2 , 
sin 2 2^i 2 = 0.7, and small mixing angle solution Am 2 2 = 10 _5 eV 2 , sin 2 2#i 2 = 0.005. It 

At72 2 

should be noted that the oscillation term S\ 2 = sin 2 1.27 — — — L cannot be neglected in 

E 

the sub-GeV experiment of atmospheric neutrino though Am 2 2 <C Am 23 , because this 
oscillation term is not so small in sub-GeV neutrino energy (E = 0.2 — 1.3GeV) with 
zenith angle 6 ~ 180°(L ~ 10000km). 

In the atmospheric neutrino experiments, the matters of the Earth have a important ef- 
fect. Matter effect is represented by a term A = 2\f2EGpN e = 7.59pE x 10~ 5 eV 2 induced 
by the Earth matter, where p is the matter density. In the Earth, p = 3.5 - 13g/cm 3 , and 



sub-GeV and multi-GeV neutrino energy ranges from O.lGeV to lOOGeV, then the value 
of A goes through those of Amf 2 and Am| 3 and the resonance happens in mixing angles 
9 12 and 9 13 as will be seen in Eqs.(8)-(11) and Eqs.(13)-(16). In the sub-GeV neutrino 
energy, A sin 26 , 13 /2Am2 3 -C 1, then we can approximate the mixing matrix U as 

U™ h = exp(iA 7 23 ) exp(zA 5 13 ) exp(zA 2 0#), (7) 

where 



sin2^ = — -^sin20 12 , 



Am? 
Am] 2 

Amff = mf 2 - mf 2 , E = m 2 + m 2 , (9) 



mf 2 2 = ^{(S + Acos 2 ^ 13 ) 

=F^/(A cos 2 13 - Amf 2 cos 20 12 ) 2 + (Am? 2 sin 26> 12 ) 2 | , (10) 
m ^ 2 = m 2 + Asin 2 13 . (11) 

On the other hand, J 4sin2# 13 /2Am 23 is not so small in multi-GeV neutrino and upward 
muon energy In this case, the mixing matrix can be expanded by powers of a quantity 
Am 2 2 sin 2^ 12 /2m| which is very small. In the lowest order, the mixing matrix is expressed 
as 

U^ ti = ex P (iA T 23 ) exp(zA 5 ^), (12) 

where 

/\ 777 2 

sin2 ^ = w^ sin2 ^ 13 ' (13) 

Amff = mf 2 - mf 2 , A = E - Am 2 2 cos20 12 , E = m\ + m 2 , (14) 

T\J (A - Am 2 3 cos 20 13 ) 2 + (Am 2 3 sin 26» 13 ) 2 } , (15) 
m f 2 = i(S + Am 2 2 cos 20i 2 ). (16) 

These expressions correspond to the neutrino case, and A for the anti-neutrino case has 
an opposite sign of the neutrino case. From these expressions, we can recognize that 



T in (M k ,x k ) = J2(U Mk )i'a(U Mk *)iaexp (- 

a \ 



a resonance exists at A m Amj 3 cos 29\3 of sin 29f£. When Amf 3 = 3 x 10 _3 eV 2 and 
cos2^ 13 ~ 1, the resonance occurs at E ~ 10 GeV. 

When the density of the Earth matter is constant, the transition probability P M [y\ — > 
vy) with the matter effects can be expressed by the expression (|5|) with U and Am^ 
replaced by U M and Amf^ 2 . Actual density of the Earth is not constant but the density 
in shells which compose the Earth is almost constant. The net transition probability 
Pvi-*v,, of neutrinos going through the Earth can be given by connecting the transition 
amplitudes of each shell. The transition amplitude in k-th shell is expressed as 

iH- (17) 

where U Mk represents the mixing matrix containing the k-th shell matter effects, m^f k 
the mass of u a in the k-th shell and x the length traveling in the k-th shell. Using this 
transition amplitude, we get the transition probability as follows: 

■■■T hh {M 2 ,x 2 )T hl {M 1 ,x 1 )\ 2 . (18) 

Calculation using this expression is very complicated and cannot be carry out analytically, 
then we calculate this numerically. 

3 NUMERICAL ANALYSES OF ATMOSPHERIC 
NEUTRINOS 

In this work, we analyze the ratio Ne-x. p {1) /Nyic(l) of experimentally observed events 
-^Exp(0 an d expected events N MC (l) without oscillation, where I represents //-like and 
e-like events. The zenith angle 9 dependent events dN-E xp (l) / d cos 9 and dNuc(l) I d cos 9 
are defined as 

d j^g =£/ e l {E l )a Vl {E Vin E l ^)F Vl ,{E Vin 9-^)P M {v ll - ^dE^d cos ^ (19) 
^ Mc( j } = [ e l (E l )cx Ul (E Ul ,E l ,tfj)F Ul (E Ul ,9-^)dE Ul dE l d C osif;, (20) 

& COS (/ J 
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where the summation Ylv, are taken over and v e . In these expressions, processes of z> M 
and v e are contained. ei(E{) is the detection efficiency of the detector for /-type charged 
lepton with energy E\. a Ul (E Vl , Ei, tp) is the differential cross section of scattering / with 
energy E\ by incident v% with energy E v = E VI , where angle ip is the scattering angle 
between the directions of incident v\ and scattered I. F U/ (E UI ,8) is the incident Vy flux 
with energy E V l produced at the atmosphere coming to the detector with zenith angle 9. 
E Ul and 9 dependences of F n {E vv 9) for multi-GeV experiment [E Vl > 1.33GeV) are given 
in Refs. [fU], P}| . These dependences including the geomagnetic effects for sub-GeV case 
(0.2GeV < E Vl < 1.33GeV) are taken from Ref. [12]. Other informations of ti{E{) and 
a Ul (E Ul , E[,ip) are given in Ref. (T^|. The upward through-going muons (thru-muons) and 
stopping muons (stop- muons) fluxes are shown in Ref. [|n|. Typical energies of that 
produce thru- and stop-muons are lOOGeV and lOGeV, respectively. Explicit calculation 
of Eqs. (19) and (20) is explained precisely in Appendix A of the second paper in Ref. ||. 

P M {y\i — > V]) is the transition probability with the matter effects expressed in Eq. (18). 
The Earth consists of 5 shells approximately |jl4 |; the density of the most outside shell 
(radius r = 5991 — 6371km) is 3.5g/cm 3 , the next outside shell (r = 5719 — 5991km) 
4g/cm 3 , the middle shell (r = 3519 — 5719km) 5g/cm 3 , the outer core shell (r = 1231 — 
3519km) llg/cm 3 , the inner core (r = — 1231km) 13g/cm 3 . A number of the shells 
through which neutrinos pass changes with a change of the zenith angle 9. For example, 
the number of the shell is for < 9 < 90° and 9 for 9 = 180°. 

The zenith angle distributions of sub-GeV, multi-GeV neutrino events and upward 
thru- and stop-muons fluxes are given by the SuperKamiokande 850-920 live days experi- 
ments || . The data are tabulated in Table I. These values are taken from the experimental 
event data and Monte-Carlo simulations which are given graphically in Ref. 0. /i-like 
events include the fully contained and partially contained events. Errors represent statis- 
tical ones only. 

Since P M (i>i> — > v{) is a function of Amf 2 , Am| 3 , #i 2 , #13 and #23, the ratio {dN Exp (l) 
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Table 1: e-like, //-like atmospheric neutrino data, and upward muons fluxes of Su- 
per Kamiokande experiments (850-920 live days data). These values are the ratios of 
experimental data and Monte-Carlo simulations which are obtained from the graphs in 
Ref. 0. //-like events include fully contained and partially contained events, and errors 
represent statistical ones only. 

Sub-GeV data 

cos 9 range e-like event ratio Ne xp /Nmc //-like event ratio Ne xp /Nmc 



-1.0 — 0.8 


1.10 ±0.08 


0.66 ±0.05 


-0.8—0.6 


1.10 ±0.08 


0.52 ±0.05 


-0.6—0.4 


1.05 ±0.08 


0.59 ±0.05 


-0.4 — 0.2 


0.98 ±0.08 


0.67 ±0.05 


-0.2-0.0 


0.99 ±0.08 


0.64 ±0.05 


0.0- 0.2 


1.10 ±0.08 


0.71 ±0.05 


0.2- 0.4 


1.09 ±0.08 


0.80 ±0.06 


0.4- 0.6 


0.90 ±0.07 


0.82 ±0.06 


0.6- 0.8 


1.15 ±0.08 


0.80 ±0.06 


0.8- 1.0 


0.94 ±0.07 


0.87 ±0.06 



Multi-GeV data 


cos 9 range 


e-like event ratio N^ p /Nmc 


//-like event ratio N-£ xp /Nmc 


-1.0—0.8 


0.83 ±0.18 




0.47 ±0.07 


-0.8—0.6 


1.19 ±0.18 




0.51 ±0.07 


-0.6—0.4 


0.97 ±0.15 




0.45 ±0.06 


-0.4 — 0.2 


1.18 ±0.16 




0.55 ±0.07 


-0.2-0.0 


1.04 ±0.14 




0.63 ±0.07 


0.0-0.2 


0.94 ±0.13 




0.81 ±0.08 


0.2-0.4 


1.33 ±0.16 




0.91 ±0.09 


0.4-0.6 


0.85 ±0.15 




0.99 ±0.10 


0.6-0.8 


1.00 ±0.17 




0.85 ±0.09 


0.8-1.0 


1.27 ±0.21 




1.00 ±0.10 



Upward through-going // data 
cos Grange // flux ratio A^e xp /^mc 



-1.0- 


-0.9 


0.84 ±0.11 


-0.9- 


-0.8 


0.90 ±0.11 


-0.8- 


-0.7 


0.72 ±0.09 


-0.7- 


-0.6 


0.95 ±0.10 


-0.6- 


-0.5 


0.86 ±0.09 


-0.5- 


-0.4 


0.76 ±0.08 


-0.4- 


-0.3 


0.82 ±0.09 


-0.3- 


-0.2 


1.09 ±0.09 


-0.2- 


-0.1 


0.96 ±0.08 


-0.1 


-0.0 


1.14 ±0.08 



Upward stopping // data 
cos 9 range // flux ratio N Exp /N MC 



-1.0- 


-0.8 


0.55 ±0.10 


-0.8- 


-0.6 


0.50 ±0.09 


-0.6- 


-0.4 


0.60 ±0.09 


-0.4- 


-0.2 


0.47 ±0.08 


-0.2 


-0.0 


0.75 ±0.10 
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/ d cos 9) /(dNuc(l)/d cos 9) of the zenith angle distributions is a function of Amf 2 , Am| 3 , 
#12, #13, #23 and 9. We analyze the atmospheric neutrino data fixing the values of param- 
eters Amf 2 and sin 2 2#i 2 determined from the solar neutrino experiments || as Am 2 2 = 
6 x 10" 5 eV 2 and sin 2# 12 = 0.7, which corresponds to the large angle solution, and 
Am 2 2 = 10~ 5 eV 2 and sin 2 2# 12 = 0.005, which corresponds to the small angle solution. 
We treat the ratios of the zenith angle distributions of the experimental events and the 
ones of Monte-Cairo simulation, (A^" Exp (/)/iV MC (/)) i , where % represents the region number 
of the bins of zenith angle 9. 

We calculate numerically the x 2 defined as 

a _ {(N Exp (l)/N MC (l))f - (iV Exp (0/iVMc(0r ta } 2 

x h Mf + M 2 ' ( } 

For the sub-GeV, multi-GeV and upward thru-muon experiments, the summation on i 
are from 1 to 10 of zenith angle range bins and for the upward stop-muon, from 1 to 5. 
The summation on / are over and e for the sub-GeV and multi-GeV experiments. cr st 
represents the statistical error and a SJ systematic one. We assumed that the value of cr sy 
is 5% of an amount of the A^mc for the sub-GeV and multi-GeV neutrino data, 10% for 
the upward stop-muon data and 20% for the upward thru-muon data. We adopted the 
rather large value of a sy for upward thru-muons, because there exists larger uncertainty 
of neutrino flux for high energy neutrinos. We estimated the values of x 2 f° r the various 
values of Am 23 , #13 and #23- 

In Fig. [I], we showed the contour plots of x 2 m the tan 2 #i3-tan 2 # 2 3 plane for various 
values of Am 23 . These plots show the x 2 of the sub-GeV neutrino plus multi-GeV neutrino 
zenith angle distributions for the large 9 12 angle solution. We did not show the plots for 
the small angle solution, because x 2 f° r the small angle solution is almost similar to that 
for large one. In these figures, broken thick, broken thin and dotted curves denote the 
regions allowed in 99%, 95% and 90% C.L., respectively. From these plots we can say 
that the values of allowed Am 23 is from 2 x 10 _3 eV 2 to 10 _2 eV 2 , #13 is < 17° and #23 is 
from 35° to 55°. 
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Figure 1: The plots of allowed regions in the tan 2 # 13 -tan 2 #23 plane determined by the sub- 
GeV neutrino plus multi-GeV neutrino zenith angle distributions of SuperKamiokande 
848 live days data. These figures correspond to the large 812 angle solution, Am 2 2 = 
6 x i(T 5 eV 2 and sin 28 = 0.7. In these figures, broken thick, broken thin and dotted 
curves denote the regions allowed in 99%, 95% and 90% C.L., respectively. 
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Figure 2: The plots of allowed regions in tan 2 6> 13 -tan 2 Q 2 z plane determined by the 
combination of the zenith angle distributions of the upward thru- and stop-muon Su- 
perKamiokande 902-923 live days data for large #12 angle solution. There is no difference 
between the large #12 angle solution and small one. In these figures, broken thick, broken 
thin and dotted curves denote the regions allowed in 99%, 95% and 90% C.L., respectively. 
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Figure 3: The plots of the allowed regions in tan 2 # 13 - tan 2 # 2 3 plane determined by the 
combination of the zenith angle distributions of the sub- and multi-GeV neutrinos and 
upward thru- and stop-muon data for large 9 12 angle solution. There is no difference 
between the large 6*12 angle solution and small one. In these figures, the broken thick, 
broken thin and dotted curves denote the regions allowed in 99%, 95% and 90% C.L., 
respectively. 
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Plots of the upward thru- plus stop-muons are shown in Fig. |2| for the large On angle 
solution. We did not show the plots for small angle solution, because there is no difference 
between the large On angle solution and small one. In these figures, sharp allowed regions 
similar to those of previous sub-GeV plus multi-GeV neutrinos data is not appeared. The 
reason is as follows: although stop-muons data show the sharp allowed regions similar 
to those of sub-GeV plus multi-GeV neutrino data, thru-muons data show the excluded 
regions near 45° of #23- The latter situation is caused from small deficit (about 0.2) from 
1 of the \i event ratio N Exp /N MC as shown in Table 1. 

In Fig. 3, we showed the contour plots of x 2 f° r the combination of the sub-GeV 
and multi-GeV neutrinos, the upward thru- and stop-muons zenith angle distributions for 
large On solution. In these figures, broken thick, broken thin and dotted curves denote 
the regions allowed in 99%, 95% and 90% C.L., respectively. 

From these plots, we can get the following results for the neutrino mass and mixing 
parameters: 

(1) As shown in Fig. 1, the allowed region for Am^ obtained from the sub-GeV and 
multi-GeV neutrino experiments is from 2 x 10~ 3 eV 2 to 10~ 2 eV 2 at 90%C.L., and the 
allowed region for #23 angle is from 35° to 55° and for 13 angle is less than 17°. The 
minimum \ 2 is obtained as 29 for 39 degrees of freedom (DOF) at Am^ = 5 x 10 _3 eV 2 , 
#13 = 10° and #23 — 45°. The minimum \ 2 i n the restriction # 13 = which corresponds to 
the two flavor mixing — v T approximation, is 30 for 39 DOF at Am| 3 = 4 x 10 _3 eV 2 . 

(2) As shown in Fig. 2, the allowed region for 13 and #23 is excluded at small #13 and the 
maximal mixing of O23 for large values of Am| 3 . This is because of the fact that the data 
in upward thru-muons has not so large deficit from the no-oscillation case as seen in the 
Table 1. 

(3) There is no significant difference between the large #12 solution and the small solution. 
This fact is not seen from these figures, but we can see from our numerical calculation 
that the difference between the large angle solution and the small one is less than 1% for 
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multi-GeV neutrinos and upward muons and less than 10 % for sub-GeV neutrinos. 

(4) The allowed region for Am| 3 obtained from the sub- and multi-GeV neutrinos, and 
upward thru- and stop-muons is given as 2 x 10~ 3 eV 2 to xlO _2 eV 2 , and for 823 as 35° 
to 55° and for # 13 as less than 13°. The minimum x 2 is obtained as 55 for 54 DOF at 
Am| 3 = 4 x 10~ 3 eV 2 , 13 = 10° and 823 = 45°. The minimum x 2 in the restriction 6> 13 = 0° 
which corresponds to the two flavor mixing — v T approximation, is 61 for 54 DOF at 
Am^ = 3 x 10~ 3 eV 2 and #23 — 45°. These results are the same as the result obtained by 
Super-Kamiokande collaboration [0]. 

(5) It is interesting that the minimum of x 2 is obtained at not 813 = but 813 = 10°, though 
the difference between x 2 f° r $13 — case and that for # 13 = 10° case is not so large. This 
result is consistent with the CHOOZ experiment 0. The CHOOZE experiment predicts 
the results: B 13 < 13° for Am| 3 = 3 x l(T 3 eV 2 . 

We will now discuss about the interesting feature that the mixing angle #13 is preferred 
to be about 10°. If this feature is real, the detected v e events in the long baseline K2K 
experiment || will be about 10 times as large as the events expected in 8 13 = case. 
This is caused from the fact that the flux produced at KEK is about 100 times as large 
as z/ e flux and the transition probability is Piy^ — > v e ) ~ sin 2 2# 13 sin 2 1.27 /S.m^L/ E, 
E ~ 1.4GeV and L = 250km. 

4 Conclusion 

We analyzed the atmospheric neutrino experimental data of Super-Kamiokande |l[] in the 
three-flavor neutrino framework with the mass hierarchy m\ ~ ui2 m 3 and obtained 
the allowed regions of parameters Am^, 813 and #23, including the Earth matter effects 
thoroughly. We studied the event ratios of the sub- and multi-GeV, and upward thru- and 
stop-muons zenith angle distributions. From these atmospheric experiments, we can get 
the allowed region of mass parameter Am^ restricted as 0.002eV 2 -0.01eV 2 , and mixing 
parameter 813 as less than 13° and #23 as 35°-55°. The value of Am| 3 at the minimum 
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X 2 = 55 for 54 DOF is obtained as 4 x 10~ 3 eV 2 at 6>i 3 = 10° and 6> 23 = 45°. The minimum 
X 2 = 61 for 54 DOF is obtained with the restriction 6> 13 = at the Am^ = 3 x 10~ 3 eV 2 
and 623 = 45°. This fact seems very interesting for us because the mixing parameter #13 
may not be 0. For mixing parameter 812, the difference between the large angle solution 
and the small one is less than 1%. If # 13 is about 10°, the detected v e events in K2K 
experiment is about 10 times as large as events expected in 6> 13 = case. 
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